From Point Defects in Graphene to Two-Dimensional Amorphous Carbon 
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While crystalline two-dimensional materials have become an experimental reality during the past 
few years, an amorphous 2-D material has not been reported before. Here, using electron irradiation 
we create an sp^-hybridized one-atom-thick flat carbon membrane with a random arrangement 
of polygons, including four-membered carbon rings. We show how the transformation occurs step- 
by-step by nucleation and growth of low-energy multi-vacancy structures constructed of rotated 
hexagons and other polygons. Our observations, along with first-principles calculations, provide 
new insights to the bonding behavior of carbon and dynamics of defects in graphene. The created 
domains possess a band gap, which may open new possibilities for engineering graphene-based 
electronic devices. 

PACS numbers: 68.37.0g, 81.05.ue, 64.70.Nd, 31.15.es 



Hexagonal rings serve as the building blocks for the 
growing number of sp^-honded low-dimensional carbon 
structures such as graphene [U [2] and carbon nan- 
otubes [3]. Non- hexagonal rings usually lead to the de- 
velopment of non-zero curvature, e.g.^ in fullerenes [4] 
and carbon nanohorns [5], where the arrangement of 
other polygons can be geometrically deduced via the 
isolated pentagon rule (IPR) [H E] and Euler's theo- 
rem [7]. Aberration corrected high-resolution transmis- 
sion electron microscopy (AC-HRTEM) has recently al- 
lowed atomic-resolution imaging of regular carbon nanos- 
tructures and identification of defects in these materi- 
als [8HT2]. Point defects, mostly vacancies, are naturally 
created by the energetic electrons of a TEM. However, 
the possibility for selectively creating topological defects 
representing agglomerations of non-hexagonal rings could 
be more desirable in the context of carbon-based elec- 
tronics. ^13H15] 

In fact, despite the recent advances, the precise mi- 
croscopic picture of the response of graphene to electron 
irradiation remains incomplete. Earlier experiments on 
curved carbon nanosystems ^Jg' ^XT^ have shown that they 
avoid under-coordinated atoms under irradiation at high 
temperatures via vacancy migration and coalescence. Re- 
cent experiments on graphene [12 reported only the de- 
velopment of holes. Theoretical studies have also pre- 
dicted the appearance of small holes or formation of 
haeckelite-like configurations [18 or dislocations [19]. 

In this Letter, we report the transformation of 
graphene into a two-dimensional random arrangement 
of polygons due to continuous exposure to the electron 
beam with an energy just above the knock-on thresh- 
old. By carefully choosing the electron energy, we se- 
lectively enhance and suppress the underlying mecha- 
nisms of defect production. A combination of experi- 
ments and density-functional theory (DFT) calculations 
allows us to show that the transformation is driven by 
two simple mechanisms: atom ejection and bond rota- 



tion. The created defects tend to have a low formation 
energy and exhibit an electronic band gap. We also dis- 
cover other unexpected configurations, such as stable car- 
bon tetragons [20] which appear upon linear arrangement 
of di- vacancies. 

Our graphene membranes were prepared by micro- 
mechanical cleavage and transfer to TEM grids [21 . 
Aberration-corrected HRTEM imaging was carried out 
in an FEI Titan 80-300, equipped with an objective-side 
image corrector. The microscope was operated at 80 keV 
and 100 keV for HRTEM imaging, and at 300 keV for ir- 
radiation. The extraction voltage of the field emission 
source was set to a reduced value of 2 kV in order to 
reduce the energy spread. For both 80 keV and 100 keV 
imaging, the spherical aberration was set to 20 A/im and 
images were obtained at Scherzer defocus (ca. —9 nm). 
At these conditions, dark contrast can be directly inter- 
preted in terms of the atomic structure. 

The DFT calculations were carried out with the VASP 
simulation package [22, ^23^ using projector augmented 
wave potentials [24] to describe core electrons, and the 
generalized gradient approximation [25 for exchange and 
correlation. Kinetic energy cutoff for the plane waves was 
500 eV, and all structures were relaxed until atomic forces 
were below 0.01 eV/A. The initial structure consisted 
of 200 C atoms, and Brillouin zone sampling scheme of 
Monkhorts-Pack [26] with up to 9 x 9 x 1 mesh was used 
to generate the /c-points. Barrier calculations were car- 
ried out using the nudged elastic band method as imple- 
mented in VASP [27]. 

We started our experiments by monitoring in situ the 
behavior of graphene under a continuous exposure to 
electron irradiation using AC-HRTEM imaging with an 
electron energy of 100 keV, i.e., just above the threshold 
for knock-on damage {Td) in sp^ -bonded carbon struc- 
tures 2H1 113 • Fig- [1^ shows a graphene structure after 
an electron dose of ^ 1-10^^ e~/nm^. Contrary to the ex- 
pectations, the structure does not predominantly consist 
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Figure 1. (a) (Color online) Amorphous two-dimensional sp^-bonded carbon membrane created by a high-dose exposure of 
graphene to 100 keV electron irradiation in an HRTEM. The blue colored area is a single-layer carbon structure. Scale bar is 
1 nm. (b) Fourier transforms (power spectra) from HRTEM images of the initial graphene configuration (top), an intermediate 
configuration (center), and of the amorphous 2D carbon (bottom). 



of holes or collapse into a 3D object. Instead, it has re- 
mained as a coherent single-layer membrane composed of 
a random patch of polygons. Holes have also formed, but 
only on a small fraction of the area. The Fourier trans- 
form of the image shows that the resulting structure is 
completely amorphous (Fig[T]3). 

In order to understand the mechanisms behind the 
transformation, we separated them by varying the elec- 
tron beam energy. To observe how a defected graphene 
sheet reacts to an electron beam when atomic ejections 
are prohibited by a low enough electron energy, we cre- 
ated initial damage in a graphene sheet by brief 300 keV 
irradiation, and then studied the generated structures 
at 80 keV. Now only under-coordinated atoms can be 
ejected (T^ for a sp^-bonded C is about 18-20 eV [281129]. 
whereas DFT calculations predict a of ~14 eV for 
a two-coordinated C). However, bond re-organization 
is possible, as activation energies for bond rotations in 
sp^-bonded carbon structures are in the range of 4- 
10 eV [30, 31 , depending on the local atomic configu- 
ration. Correspondingly, in pristine graphene, bond ro- 
tations are occasionally observed under 80 keV irradi- 
ation [9^ (Fig. [2|L,b), resulting in the formation of the 



Stone- Wales defect [32 . In all observed cases continued 
exposure reversed this transition in pristine graphene. 
However, defect structures, e.g.^ di- vacancies, can con- 
vert between different configurations (Fig.[2]i-f) via bond 
rotations. According to our DFT simulations, the bar- 
riers for bond rotations in these structures are 5-6 eV, 
which excludes thermally activated migration. 

In Fig. |3^-d we present evolution of a more complex 
defect structure. The defects created by an electron 
beam are predominantly mo no- vacancies, which quickly 
convert to di-vacancies due to a higher probability for 
under-coordinated atoms to be ejected, as noted above. 
Here (Fig. [3^), a brief exposure to a 300 keV beam (dose 
^ 10^ e~/nm^) has created an isolated V2 (555-777) and 
a defect with 4 missing carbon atoms (two connected 
di-vacancies). During the image sequence, the 80 keV 
electron beam causes the structure to re-organize via 
bond rotations. The (555-777) turns first into a ^2(5- 
8-5) (Fig. [s}) and then a dislocation dipole (Fig. [s]^), be- 
fore forming a defect composed of clustered di-vacancies 
(Fig. [sji). Fig.|3^-d also present two frequently observed 
linear arrangements of di-vacancies. 

Because atom ejection occurs at random positions, va- 




Figure 2. (Color online) Elementary defects and frequently observed defect transformations under irradiation. Atomic bonds 
are superimposed on the defected areas in the bottom row. Creation of the defects can be explained by atom ejection and 
re-organization of bonds via bond rotation, (a) Stone- Wales defect, (b) defect-free graphene, (c) Vi{5-9) single vacancy, (d) 
V2 (5-8-5) di- vacancy, (e) V2 (555-777) di- vacancy, (f) V2 (5555-6-7777) di- vacancy. Scale bar is Inm. 



cancies initially appear randomly in the area exposed to 
the electron beam. However, during lower-energy expo- 
sure {e.g., 80 keV), these defects travel via a re-bonding 
mechanism, which is illustrated in Fig. [3^. Each migra- 
tion step is initiated by a single electron impact from 
which the atom obtains energy slightly below T^. In other 
words, electron irradiation provides the activation energy 
to drive the system from a local energy minimum into 
another one, in our case predominantly via (reversible) 
bond rotations (Fig.[2]i-f, Fig.[3^-d). This can be clearly 
seen in video S4 in Ref. 33 (partially shown in Fig.[3|i-d), 
where the configuration changes frequently until it ar- 
rives in the more stable configuration composed of three 
aligned double- vacancies. 

Individual transitions can also lead to higher energy 
structures. For example, the intermediate states for the 
di-vacancy migration [dislocation dipole (Fig. ^) and 
V2 (555-777) (Fig. [3^)] have formation energies Ef which 
differ from that of the I/2 (5-8-5) by - +3.72 eV and 
~ —0.66 eV, respectively. Aligning di-vacancies along 
the zigzag direction of the lattice (Fig. [3]i), ~ 1.32 eV 
of energy is gained per a di-vacancy pair, as compared 
to isolated di-vacancies. When the alignment appears 
along the armchair direction (Fig.|3^-c), the energy gain 
is 2.01 eV. In this case, a tetragon is formed where the 
two pentagons of the adjacent di-vacancies would overlap. 
HRTEM simulation of the DFT-optimized structure of 
the defect is in excellent agreement with the experimen- 
tal image [33 . Note that sp^-bonded carbon tetragons in 
molecules, as in cyclobutadiene, can be stabilized only at 
low temperatures and when the molecules are embedded 
into a matrix [20 . In our case they are stabilized by the 
surrounding graphene lattice, as theoretically predicted 
for nanotubes [34 . 

Under 100 keV irradiation, atom ejection occurs at a 
very slow rate, so that changes in the atomic network 
are sufficiently slow to be precisely resolved. Therefore, 



the reconstruction of vacancy defects via bond rotations 
can be monitored immediately after a vacancy is gen- 
erated. An example is presented in Fig. [4j The initial 
configuration (Fig. |4^) consists of three di-vacancies in 
the armchair orientation (formed prior to recording the 
first image). In the recorded images, the structure loses 
atoms until 24 atoms are missing. Several remarkable 
configurational changes are found in this image sequence. 
Fig. |4^-c shows a collapse of linearly clustered defects 
into an apparently less defective structure with a dis- 
location dipole. This corresponds to the prediction of 
Jeong et al. [19 that the dislocation dipole is favored 
over a large mult i- vacancy. However, this requires a lin- 
ear arrangement of vacancies. Fig. ^-d shows the loss of 
four additional atoms. Two of them gave rise to an addi- 
tional di-vacancy; the other two contributed to the sepa- 
ration of the dislocation cores (the rotated hexagon, clus- 
tered with a Stone- Wales defect, constitutes a dislocation 
core). During continued irradiation, we see the formation 
of a cluster of rotated hexagons surrounded by a chain of 
alternating pentagons and heptagons (Fig.[4^-h). Such 
configurations, rotated by 30° with respect to the origi- 
nal lattice and matched by pentagons and heptagons to 
the zigzag lattice direction, appear to be the preferred 
way to incorporate missing atoms in the graphene struc- 
ture. Due to the matching numbers of pentagons and 
heptagons - and hence cancellation of negative and pos- 
itive curvature - these structures remain flat. 

To understand the driving force for the transforma- 
tions, we calculated Ef for the simplest defect struc- 
tures matching the obverved trend of forming a rotated 
hexagon kernel. The lowest-energy tetra-vacancy (four 
missing atoms) can be created by combining two ^2(5555- 
6-7777) di-vacancies, whereas the hexa- vacancy (six miss- 
ing atom) requires three of these defects (Fig. [5|i). Re- 
markably, these configurations have the lowest Ef of any 
reported vacancy structures with equal number of miss- 
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Figure 3. (Color online) (a-d) Electron beam driven di- vacancy migration observed at 80 keV. In (e), the changes in the bond 
configuration required for allowing the migration of a di- vacancy are shown. Transformation 1/2 (5-8-5) ^2 x (5-7) is initiated 
by rotating bond A, and V2(5-8-5)^V2 (555-777) by rotating bond B. Rotating bonds C and D will lead to the final 1^2(5-8-5) 
structures. In the first case the defect has moved by ai — a2, and in the second case it has rotated by 60° around pentagon 2. 
The original TEM images without overlays for panels (a-d) are presented in Ref. [33] . 



ing atoms in graphene {Ef per missing atom multiplied 
by the number of missing atoms are ~ 4x3.14 eV and 
~ 6x2.50 eV). A hole with six missing atoms has a for- 
mation energy of ^ 6x3.15 eV, while for a dislocation 
and a haeckelite-like structure values of ~ 6x3.72 eV and 
~ 6x2.64 eV have been reported [T^, respectively. Ev- 
idently, the rotated hexagon defects spawn the family 
of lowest energy mult i- vacancies in graphene. In con- 
trast to what was recently shown for the zigzag-oriented 
di- vacancies [14 , these structures open a band gap in 
graphene, as can be seen from the density of states 
(Fig. [sJd). The calculated band gap is in the order of 
200 meV. This value is possibly underestimated within 
the used GGA approximation, and advanced DFT meth- 
ods are likely to give a higher value j35] . 

To conclude, we have shown how an electron beam can 
be used to selectively suppress and enhance bond rota- 
tions and atom removal in graphene. We demonstrated 
that irradiation at electron energies just above the 
threshold for atom displacement turns graphene not into 
a "perforated graphene" but a two-dimensional coher- 
ent amorphous membrane composed of sp^-hybridized 
carbon atoms. This membrane grows through nucle- 
ation and expansion of defects which result from elec- 
tron beam-driven di-vacancy migration and agglomer- 
ation. These defect configurations predominantly con- 
sist of a 30° rotated kernel of hexagons surrounded by a 
chain of alternating pentagons and heptagons. These de- 
fects are the energetically favored way for the graphene 



lattice to accommodate missing atoms, and have a semi- 
conducting nature. Since several of the presented exam- 
ples of the two-dimensional sp^-hybridized defect con- 
figurations violate the IPR, due to increased reactiv- 
ity [6l |36], they may be exploited for functionalization 
of graphene. We also showed unambiguous evidence 
for four-membered carbon rings in graphitic structures. 
Clearly, despite of the large amount of research, the rich- 
ness of carbon chemistry continues to provide surprises. 
More examples of the observed structures and videos of 
the complete TEM image series are presented in Ref. [33] . 
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Figure 4. (Color online) Formation of rotated hexagon-kernels in mult i- vacancy structures under a 100 keV electron beam. 
Scale bar is Inm. The original TEM images without overlays are presented in Ref. [33] . 




Figure 5. (Color online) Idealized rotated-hexagon defects (a) formed from one, two and three V2 (5555-6-7777) di- vacancy 
defects, as optimized with DFT calculations. Electronic density of states (b) for pristine graphene and the structures presented 
in panel (a). 
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METHODS 

As described in the main article, graphene membranes 
were prepared by micro-mechanical cleavage and transfer 
to TEM grids [21 . Aberration-corrected HRTEM imag- 
ing was carried out in an FEI Titan 80-300, equipped 
with an objective-side image corrector. The microscope 
was operated at 80 keV and 100 keV for HRTEM imag- 
ing, and at 300 keV for irradiation. The extraction volt- 
age of the field emission source was set to a reduced value 
of 2 kV in order to reduce the energy spread. For both 
80 keV and 100 keV imaging, the spherical aberration 
was set to 20 fim and images were obtained at Scherzer 
defocus (ca. —9 nm). At these conditions, dark contrast 
can be directly interpreted in terms of the atomic struc- 
ture. Image sequences were recorded on the CCD camera 
with exposure times ranging from 1 s to 3 s and intervals 
between 4 s and 8 s, and a pixel size of 0.2 A. The effect of 
slightly uneven illumination is removed by normalization 
(division) of the image to a strongly blurred copy of the 
same image, effectively removing long-range variations. 
Drift-compensation is done using the Stackreg plugin for 
the Image J software [37]. We show individual exposures 
as well as averages of a few frames (up to 10). This is 
because we have used different beam current densities, 
in order to test for possible dose rate effects (within the 
100 keV, simultaneous imaging and defect generation ex- 
periment). The configurations described here can be dis- 
cerned in individual exposures if ca. 10^ counts per pixel 
at 0.2 A pixel size are used. With a high beam current 
density, this is possible in 1 s exposures (corresponding 
to a total dose of ca. 10^ per image, and a dose rate 

of 10^ ^4—). At lower current densities, the same dose 
was spread over several exposures, so that sample drift 
could be compensated (individual exposures could not be 
longer than 3 s due to sample drift). The lowest dose rate 
was ca. 3 • 10^ ^|—. Within our dose rate range, density 
and shape of defects appears to depend only on the total 
dose. 

Formation energy of a structure in our DFT calcula- 
tions was defined in the usual way as 

Ef{V„) = E,ot{Vn) - ^fot + nfi^r, (1) 

where £^tot(K,) and Ef^^ are the total energies of the 
structure with the defect (n missing atoms) and the same 
supercell without the defect, respectively, and figr is the 
chemical potential of a carbon atom in pristine graphene. 
The semi-conducting features of the defects were consis- 
tently observed with varying number of /c-kpoints used 
in calculating the electronic density of states. 



FURTHER IMAGES 

For clarity we show here the TEM images of the main 
article without stucture overlay, and also a few addi- 
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tional TEM images. Fig. [6] shows several di- vacancy 
(DV) defects in linear alignment. In particular, the car- 
bon tetragon is always reproduced when the DVs are 
aligned along the armchair direction of the graphene lat- 
tice. Fig. [7| shows a comparison between a HRTEM im- 
age simulation based on a DFT-optimized structure and 
an actual HRTEM image for the armchair alignment of 
di- vacancies, which contains the carbon tetragons. 

In Fig. [8) we show the panels a-d from the Figure 3 of 
the main article with and without overlays. These images 
are also contained in Supplementary video S2. In Fig. |9] 
we show panels a-h from the Figure 4 of the main article 
with and without overlay (the corresponding time series 
is shown in Supplementary video S3). 

In Fig. [lOj we show two additional images where a de- 
fect with a rotated hexagon kernel was generated from 
clusters of multiple vacancies. These configurations were 
created by a short 300 keV exposure and subsequent 
imaging at 80 keV. It should be noted that the rotated 
hexagon kernels appear in the larger vacancy clusters in 
both of our experiments (100 keV irradiation with simul- 
taneous imaging as in Figure 4 of the main article, and 



the 300 keV/80 keV combination as shown in Fig. 10). 



DESCRIPTION OF THE SUPPLEMENTARY 
VIDEOS 

Supplementary video SI: Generation and transformation 
of defects under 100 keV irradiation. Each frame is 
an average of 2 CCD exposures. Note the continuous 
increase in defect density and partial amorphization of 
the membrane. 



Supplementary video S2: Further transformation of 
graphene to a 2D amorphous membrane under 100 keV 
irradiation. In terms of total dose, it can be considered 
as a continuation of video S4 (although sample region is 
different). Each frame shows an individual CCD expo- 
sure. Note that the entire clean graphene area becomes 
amorphous while remaining one-atom thick, and beam- 
generated holes make up only a small fraction of the area. 

Supplementary video S3: Two isolated di-vacancies, 
generated by brief 300 keV irradiation and observed 
at 80 keV. The di-vacancies migrate under the beam 
and transform between the (5-8-5), (555-777) and 
1^2(5555-6-7777) configurations multiple times. The 
video shows individual exposures in each frame. 

Supplementary video S4: A cluster of several vacancy 
defects, generated by brief 300 keV irradiation and ob- 
served at 80 keV. The video shows individual exposures 
in each frame. The configuration changes continuously 
until the final, linear aligned di-vancy configuration is ob- 
served and stable throughout several exposures (frames 
22-29). Frames 30-36 are duplicates of frame 29 in order 
to show the final configuration as a still image at the end. 

Supplementary video S5: Generation and transformation 
of defects under 100 keV irradiation. Each frame is an 
average of 10 CCD exposures (recorded at a lower current 
density). This video shows the generation of the rotated 
hexagon kernel as in Figure 4 of the main article. The 
rotated hexagon kernel is highlighted as overlay in frames 
28-31. 



Figure 6. Images of di- vacancy (DV) defects in linear alignment, shown with and without structure overlay in each case, (a) 
Three DVs aligned along the zigzag direction of the graphene lattice, (b) Five DVs aligned in zigzag direction, clustered with 
a single vacancy (black), (c) Two DVs aligned along the armchair direction, forming a carbon tetragon at their intersection, 
(d) Three DVs aligned along the armchair direction, forming two carbon tetragons. All scale bars are 1 nm. 




Figure 7. DFT-optimized structure of the defect with two tetragons (a), HRTEM simulation based on the DFT-structure (b) 
and the HRTEM image of the same structure for comparison (c). Scale bar is 1 nm. 




Figure 8. Panels a-d from Figure 3 of the main article with and without overlays. Scale bar is 1 nm. 




Figure 9. Panels a-h from Figure 4 of the main article with and without overlays. Scale bar is 1 nm. 
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Figure 10. Two examples of defect structures created by short 300 keV exposure and imaged at 80 keV. In the upper row 
(panels a,b) DFT-optimized configurations are superimposed over the experimental images. Lower row (panels c,d) displays the 
structures without DFT-overlays. The first structure (a,c) can be constructed by combining three V2 (5555-6-7777) reconstructed 
di- vacancies (6 missing atoms), whereas the second one (b,d) can be formed from two V2 (5555-6-7777) 's, two V2 (5-8-5) 's, one 
V2 (555-777) and one (55-77) defect (total of 10 missing atoms), as indicated by the colors. Scale bars are 1 nm. 



